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Q^ : 1. INTRODUCTION 

O; 

0^ ! Collective flow is a motion characterized by space-momentum correlations of dynamic 

Q I origin. It is of interest in collisions because it may tell us about pressures generating 

O ■ that motion and about the equation of state (EOS) and other properties of the strongly- 

^ ■ interacting matter. The flows that have been identified thus far are radial, sideward, and 

I . elliptic. 

CJ I The role played by the pressure p in the dynamics may be inferred by examining the 

rH I relativistic Euler equation in the frame of the matter, defined by vanishing collective 

L| ■ velocity v = 0: 
• 1—1 , 

rN . (p -L -n^ 



S 



e + p) — v = -Vp. (1) 

Here e is energy density. This equation is an analog of the Newton equation and enthalpy 
w = e + p plays the role of mass. Pressure plays the role of potential for the collective 
motion. 

Recent interest in the flows and EOS is due to the possibility of identifying the phase 
transition to quark-gluon plasma. Figure |1| shows a schematic diagram of temperatures 
T and baryon chemical potentials /i explored in central heavy-ion reactions with the 
indicated possible crossings of the transition at different accelerators. 

In terms of baryon density p, the pressure is: 



p = p^ ^^''P^ 



dp 



1 

s/p 



(2) 



where s is the entropy density. Expected changes in the energy per baryon e/p, if there 
were a phase transition at T = 0, are qualitatively indicated in Fig. |^. The arrow in the 
figure represents the Maxwell construction. 

Within and above any first-order phase-transition, or a transitional region, a physical 
system exhibits a softening, i.e. the ratios p/e and p/w drop. This is due to the fact that 

*e-inail: danielewicz@nscl.insu.edu 




normal matter ---> 923 1200 !^ [GeV] 



Figure 1 . Temperatures T and baryon chemical potentials /i explored in central heavy-ion 
reactions. 
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Figure 2. Energy per baryon e/p as a function of baryon density p, in units of normal 
density po, for a soft EOS (S), hard EOS (H), and a soft EOS with a high-density phase 
transition (S w/PT). 



the pressure p changes continuously across the transition, while e jumps together with 
entropy s and density p. For example, in the baryonless matter we have: 

dp = s dT , de = T ds . (3) 

With an increase in temperature across a transition, the entropy jumps and so does the 
energy. However, for the pressure only the derivative changes. Analogous changes across 
a T = transition are deduced in Fig. ^ Following (|l|), the signature of a phase transition 
should be a slowed down pace of the development of the nuclear collective motion within 
and above the phase transition. 

Equilibrium, to which EOS and our considerations so far refer to, is never precisely 
reached in reactions. This creates both difficulties and opportunites. On one hand, 
extrapolations must be done to the equilibrium limit. On the other hand, the flow may be 
exploited in exploring the transport properties of hadronic matter, related to in-medium 
cross sections. The flow may be studied to a variable degree of detail. A study of nucleon 
flow as a function of transverse momentum, for example, may reveal the dependence of 
the mean field (MF) felt by nucleons on momentum [Q. 

We now turn to the different types of flow. 

2. RADIAL EXPANSION 

The collective radial expansion is often assessed by looking for deviations of momentum 
distributions, especially transverse, from thermal. The momentum distributions are com- 
monly described in terms of the simple Siemens-Rasmussen [0] formula, or its derivatives, 

e ^ oc e--/^ ( ^ (76 + T) sinh ^ - cosh^ 1 . (4) 

dp {pv T T ) 

A safer assessment of the radial flow is by comparing spectra, or average energies, of 
particles with different mass. The higher the mass, the stronger is the effect of collective 
expansion and flatter the distribution (more spread-out by the collective velocity), and 
higher the average energy. Transverse distributions flattening with the particle mass 
have been seen in central collisions of heavy nuclei at beam energies ranging from below 
100 MeV/nucleon to above 100 GeV/nucleon, cf. Fig. ^ The larger-mass distributions are 
sharper in the lighter than in the heavier system, at the high energy, indicating a weaker 
collective expansion in the lighter system. 

Figure ^ displays excitation function of transverse temperature and of velocity in heav- 
iest systems. The velocity saturates at AGS energies, possibly due to meson production 
and progressing transparency. 

Of interest is the possible use of the radial expansion in the determination of EOS. It 
must be remembered that the separation into the collective and thermal energies occurs 
at freeze-out, when collisions become infrequent. Let us consider flrst the situation at low 
energies. If the EOS is soft and pressure low, the expansion is slower than for a stiff EOS, 
but then one just needs to wait longer for same observable values to emerge at freeze-out. 
To tell the difference, one needs some timing device. As such devices might serve the 
persistence of longitudinal motion at high energies or the early strange particle emission. 
Till now though, their exploitation in the EOS determination has been limited P]. 
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Figure 3. Transverse momentuin distributions in central collisions at 250 MeV/nucleon 
(left) and at 158 GeV/nucleon (right) from the measurements of the FOPI and of the 
NA44 [Q Collaborations, respectively. 



3. SIDEWARD FLOW 

Sideward flow is a deflection of forwards and backwards moving particles, away from 
the beam axis, within the reaction plane. The situation in reactions is schematically 
illustrated in Fig. |^. For the compressed and excited matter in a central region it is easier 
to get out to the vacuum on one side of the beam axis than on the other. Eagerness to get 
out will be enhanced by high generated pressure but also by the momentum dependence 
of MFs before the equilibration takes place. The ability to get out depends on the inter- 
particle cross sections. 

The sideward flow is often represented in terms of the mean in-plane component of 
transverse momentum at a given rapidity, {p^{y)), and additionally quantified in terms of 
the derivative at the midrapidity, see Fig. 13: 



dy 



or 



d{p^) 

diy/ys) 



(5) 



The normalization of the rapidity to the beam in the derivative enhances, somewhat 
artificially, the strength of dynamic effects at high energies relative to low. 

In transport models, it is directly observed that the production of sideward flow is 
shifted towards the high density phase as compared to the radial flow p . The sideward 
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Figure 4. Excitation functions of transverse temperature (upper panel) and transverse 
collective velocity (lower panel) at midrapidity in heavy systems. 



flow thus has more potential in the EOS determination than the radial flow. The flow 
excitation function is represented in Fig. ^ and the flow is seen to be maximal between 
0.1 and 10 GeV/nucleon. 



4. SECOND-ORDER OR ELLIPTIC FLOW 

The elliptic flow is typically studied at midrapidity and quantifled in terms of f 2: 



V2 = (cos 20) 



Vn = (cosn0) 



(6) 



where is the azimuthal angle relative to the reaction plane. The second-order flow 
may offer a better chance for the EOS determination than the flrst-order sideward flow, 
because it involves less of the uncertainties in the opposing streams of matter moving past 
each other. Typical azimuthal patterns at midrapidity may be seen in Fig. ^. 

At AGS energies the elliptic flow results from a competition between the early squeeze- 
out when compressed matter tries to move out in the direction perpendicular to the 
reaction plane and the late-stage in-plane emission associated with the shape of the par- 
ticipant zone ||lO|, cf. Fig. |[ The squeeze-out contribution to the elliptic flow depends, 
generally, on the pressure p built-up early on, compared to the energy density e, see (0), 
and on the passage time for the spectators. When the heated matter is exposed to the 
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Figure 5. In-plane particle Figure 6. Average in-plane transverse momentum compo- 
deflection. nent as a function of rapidity in central Au + Au collisions 

at 10 GeV/nucleon 0. 
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Figure 7. Excitation function of sideward flow in central collisions of heavy nuclei. 
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Figure 8. Aziinuthal distributions, with respect to the reconstructed reaction plane, 
of protons emitted from semicentral Au + Au colhsions in the rapidity intervals of 

-0.7 < y/VBeam, < -0.5 (a), -0.5 < y/VBeam < -0.3 (b), -0.1 < y/VBeam < 0.1 (c), 

0.3 < y/yBeam < 0.5 (d), and 0.5 < y/yBeam < 0.7 (e). 




Figure 9. Collision of two Au nuclei at relativistic energies. Time shots are shown for 
an instant before the collision (a), early in the collision (b), and late in the collision (c). 



vacuum in the transverse direction, a rarefaction wave moves in putting the matter into 
motion. The time for developing the expansion is then R/cg, where R is the nuclear 
radius and Cg = ^y{dp/de)s/p is the speed of sound. The passage time for spectators, 
on the other hand, is of the order of 2R/{'JqVo), where vq is the spectator cm. velocity. 
The squeeze-out contribution to the elliptic flow should then reflect the time ratio 

(7) 



70-^0 



The result (|^) gives hope that signiflcant changes in the dependence of pressure on 
energy density, such as associated with a phase transition,^ might be revealed in the 

^Discussion of the physics of refraction into mixed phase or across the transition may be found in [[ll| . 



variation of the elliptic flow. Overall, the squeeze-out contribution should decrease as 
a function of energy with the flow becoming positive, f 2 > 0. The elliptic-flow excitation 



function is displayed in Fig. 10. It is seen that, indeed, while dominated by squeeze-out 
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Figure 10. Elliptic- flow excitation function. 



at moderate energies, the elliptic flow becomes positive at high energies. Whether or not 
any changes in V2 with energy might be associated with any phase transition requires 
comparisons to transport-model calculations to assess the magnitude of such possible 
changes. It should be noted, due to the expansion of elliptically shaped participant zone 
exposed from all sides to the vacuum. Fig. ^(c), the in-plane elliptic flow is expected to 
be present at RHIC energies, see Fig. ^. 

To assess the practical quantitative level of sensitivity of flow observables to EOS and 
to other characteristics of the medium, we now turn to transport model simulations. 



5. TRANSPORT-MODEL COMPARISONS 

We shall examine results of a model formulated within relativistic Landau theory, with 
nucleon, pion, delta and A^* degrees of freedom. The Wigner distribution functions / of 
those particles follow a set of Boltzmann equations, of the same general for relativistically, 
as nonrelativistically. 



dt dp dr 



dr dp 



(8) 



The single particle energies e are functional derivatives of an energy functional specified 
in terms of the Wigner functions /. Specific parametrizations of the functional lead to 
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Figure 11. Elliptic-flow coeflicient for nucleons and pions from RQMD simulations of 
6 = (5 - 10) fm Au + Au collisions at v^i = 200 GeV pj. 



different EOS, such as in Fig. H, and different optical potentials 

U°P\p) =e{p)-m-T{p), 



(9) 



where T is kinetic energy, that either depend or do not on the momentum p in the medium 
frame. 

Typical single-particle observables from the model [|I^], for different EOS, such as trans- 
verse momentum or rapidity distributions, compare favorably to data [1^. Figure |1^, with 
baryon density contour plots from a simulation of the 2 GeV/nucleon Au + Au reaction 
at 6 = 8 fm, helps to understand the generation of elliptic flow. The simulation was done 
assuming a hard EOS and a momentum dependence in the baryonic MFs. The partic- 
ipant matter gets compressed in Fig. |T2] to < 3 p/po and then explodes when spectator 
pieces are still present. When the simulation is done in the cascade mode of the model, 
the participant matter gets compressed to higher densities and it decompresses in quite 
a sluggish manner. At the completion of the decompression the spectator pieces are far 
gone. Figure 13 shows next the evolution of the midrapidity elliptic anisotropy in the 
two simulations. In either of the simulations, the anisotropy begins to grow in the in- 
plane direction. In the hard-EOS calculation, the growth is rapid and then the anisotropy 
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Figure 12. Contour plots of baryon density in the 2 GeV/nucleon Au + Au reaction 
at 6 = 8 fm in the reaction plane (bottom panels) and in the plane through the system 
center orthogonal to the beam axis (top panels). The simulation was done assuming a 
hard EOS and momentum-dependent U°^^. The density contours are shown at p/po = 0.1, 
0.5, 0.8 and then at 0.4 intervals. Every third contour, i.e. for p/po = 0.1, 1.2, and 2.4, 
is represented by a dashed line. The remaining contours are represented by solid lines. 
Numbers in the figure indicate time in fm/c. 



changes its sign. In the cascade calculation the growth is slow and it continues till late 
times. These findings may be understood as follows. The expansion of the participant 
region first grows primarily in the in-plane direction as in that direction the pressure gra- 
dients are the largest |T0[ . If the expansion is rapid, though, the matter encounters the 
spectator pieces on its way. Then the in-plane expansion stalls, while the expansion in the 
transverse direction can develop. On the other hand, the expansion is slow, the spectator 
pieces do not present much inhibition and then expansion can continue primarily in the 
in-plane direction untill the complete decompression. 

The presence of the spectator pieces allows then to assess the pace of expansion of the 
participant matter from the elliptic flow for participants. Figure IT^ shows the elliptic 
flow excitation functions for different calculations. The flow turns more in-plane as the 
beam energy increases and the spectator pieces move faster away from the reaction zone 
and get shortened by the Lorentz contraction. For the cascade calculation, for which 
the expansion is sluggish, the flow is of a primarily in-plane character at the displayed 
energies. The harder the EOS, the more out-of-plane the flow gets. However, there is more 
sensitivity to the EOS at the low energies when the spectator pieces remain still close to 
the participant zone than at high. In addition to the EOS, the momentum dependence of 
MFs felt by the emerging baryons affects their flow. The momentum dependence makes 
the flelds felt by high-momentum particles more repulsive, speeding up the expulsion of 
these particles from the central participant zone. 

Next question is whether a phase transition could sufficiently change the expansion as 
a function of energy to yield an observable trace in the excitation function for ellipticity. 
The answer is provided in Fig. |15| that shows the excitation function for a hard EOS with 
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Figure 13. Temporal evolution of 
the midrapidity elliptic anisotropy in 
2 GeV/nucleon Au + Au collision at 
6 = 8 fm in the cascade mode of the 
model (open circles) and for a hard EOS 
with momentum dependence in the MFs 
(filled circles). 



Figure 14. Calculated elliptic flow excita- 
tion functions for Au + Au reactions. Pan- 
els (a) and (b) show, respectively, the functions 
obtained without (NMD) and with (MD) the 
momentum dependent forces. The filled cir- 
cles, filled diamonds, and open squares indi- 
cate, respectively, results obtained using a soft 
EOS, a stiff EOS, and by neglecting the MF. 
The straight lines show logarithmic fits. 



a phase transition which softens that EOS at higher densities similarly to the case in 
Fig. Ij. At lower energies the excitation function in Fig. |15| follows the hard EOS function, 
but then it switches to the soft EOS function as the densities reached in the reactions 
increase. 

Now the question is what kind of expansion scenario is followed in Nature. Comparing 



the data in Fig. 10 to the cascade calculations in Fig. 14 at 1-4 GeV/nucleon, it is apparent 
that the expansion is much more rapid than in the cascade model at the reached densities. 
Figure |16| shows next a collection of the data with momentum-dependent calculations 
superimposed. It is seen that at 2 GeV/nucleon the expansion is as rapid as for the 
hard EOS with the momentum dependence, represented in Fig. |12l At higher energies. 



though, the expansion becomes more such as for the soft EOS with behavior as a function 



of energy strongly reminding that in Fig. |I5] for the phase transition. 

Below 1 GeV/nucleon different flow data favor a soft EOS with momentum dependence, 
e.g. |T^. Sideward data favor that EOS in fact up to 10 GeV/nucleon; admittedly, though, 
the sideward flow is more exposed to the uncertainties in stopping than the elliptic flow. 
Combining older and the newer higher-energy elliptic and sideward flow results, we can 
draw broad boundaries for the dependence of pressure on baryon density (or on energy 



density) as exibited in Fig. |T^. The pressure is certainly higher than in the cascade model; 
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Figure 15. Calculated elliptic flow excita- 
tion functions for Au + Au. The diamonds 
represent results obtained with a stiff EOS. 
The circles represent results obtained with 
a stiff EOS and with a second-order phase 
transition. The lines guide the eye. 



Figure 16. Elliptic-flow excitation func- 
tion for Au + Au. The data are from 
Refs. 0, [0, and |TB|. Calculations have 
been carried out for the soft and the stiff 
EOS with momentum dependence. 
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Figure 17. Shadowed region shows boundaries on pressure vs baryon density in T = 
symmetric nuclear matter, from flow observables. 
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Figure 18. Proton flow magnitude as a function of beam energy; the lower right panel 
shows the measured Fy, while the other three panels show identical measurements of the 
parameter F, with different transport model calculations superimposed. The error bars 
include systematic uncertainties. 



other than in the 1-2 GeV/nucleon region (vicinity of the upper boundary kink in Fig. |T^ 
there is no indication of any preference for the hard EOS. A set of sideward flow data with 
model predictions superimposed is shown in Fig. 0. Contrary to Fig. |I6|, no evidence for 
a sudden softening of EOS is seen here. As to conclusions on pressure vs baryon density, 
those drawn for asymmetric matter from neutron star properties generally pertain |TB| to 
lower densities than in Fig. 0. 

Deducing the properties of equilibrium matter requires restricting the range of the de- 
pendence of baryonic MFs on momentum. This can be done by examining the dependence 
of elliptic flow on momentum at the energies of interest. 

6. DIFFERENTIAL STUDIES OF FLOW 

Studies of flow in narrow intervals of transverse momentum, rapidity and impact pa- 
rameters can help to disentagle the effects of EOS and other characteristics of strongly 
interacting matter. Figure |1^ shows the dependence of V2 on impact parameter in Au 
-|- Au collisions at 400 MeV/nucleon. At lower impact parameters, the effects of the 
stiffness of EOS and of mean-field momentum dependence compete with each other. 
However, at higher impact parameters the predictions for the momentum-independent 
and momentum- dependent MFs well separate. Additional sensitivity to the momentum 



dependence might be expected while concentrating on flow at high momenta [^ 
Figure |2^ compares the ratio of out-of-plane to in-plane proton yields 



14 



o 

II 

* 



0.20 



0.15 



0.10 



0.05 



0.00 







400 MeV/nucleon Au + Au 


- 


• HM 


- 


" o 0M '* " 


- -H A 




- 


* 8 . 


- 


— '0 


— 


i ^ 


_ 


S S m 


_ 


; +i° ■•* 


- 


- +-0%' » i 


— 


: 8 ° 


] 


-,I..i.,,,i,,..i,,,,i...,' 


1 



400 MeV/nucleon Bi + Bi 



0.0 2.5 5.0 7.5 10.0 12.5 

b [fm] 



b = 9fm 



m*/m=0.65 



2 — 




200 400 600 800 1000 

p^ [MeV/c] 



Figure 19. Negative of the ellipticity co- 
efficient at midrapidity, as a function of 
the impact parameter, from simulations of 
Au + Au colhsions at 400 MeV/nucleon. 
The squares and circles represent, respec- 
tively, the results for MFs without and 
with momentum dependence. The ffiled 
symbols are for the incompressibility K = 
380 MeV and the open symbols are for K = 
210 MeV. 



Figure 20. Out-of-plane to in-plane ratio 
Rn = (1 — 2 f 2)/(l + 2 V2) for protons emit- 
ted at midrapidity from 400 MeV/nucleon 
Bi + Bi reactions, as a function of trans- 
verse momentum. Symbols represent data 



of Ref. |T9| while lines represent transport- 
model calculations for different momentum 
dependencies in the MF. 



R 



N 



N{90°) + N{270°) _ 1 - 2f2 
A^(0°) + A^(180°) ~ l + 2f2 



(10) 



from KaoS measurements ||T9| to the results of simulations. It is apparent that the data 
permit tuning the MF momentum dependence as labeled by the effective mass. 



Figure ^ compares next the transport-model parametrizations of MF to the microscopic 
Dirac-Brueckner-Hartree-Fock (DBHF) calculations with Bonn-A interaction. Within the 
tested region of densities there is a good agreement between the parametrizations that best 
describe data []T^ and the microscopic calculations. Other microscopic theories exhibit 
a variable degree of agreement with the optimal parametrizations. 



7. CONCLUSIONS 

Flow observables are important tools for investigating properties of hadronic matter 
in energetic collisions. Measurements of second order and first order flow at beam en- 
ergies > 2 GeV/nucleon show that the matter is significantly stiffer than in the cascade 
model. Lower and also upper bounds on the pressure may be deduced as a function of 
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Figure 21. Optical potential in nuclear matter as a function of nucleon momentum, at 
different densities, from the DBHF calculations with the Bonn-A interaction [^ (short- 
dashed and short-dash-dotted lines) and in the transport-model parametrizations (solid 
and long-dashed lines). The thin horizontal dashed lines indicate the zero value for the 
potential. 



density. Dependence of the elliptic flow on momentum probes the dependence of mean 
field on momentum at supranormal densities and may be used to constrain microscopic 
theories. Some level of softening may occur in the EOS at the energies between 2 and 
4 GeV/nucleon, corresponding to a compression of 3-4 po- Further insights could be 
gained by studying the 6-dependence of V2. 
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